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On the Currents induced in a Conductor by the Passage of a 

Mass of Magnetic Material over it. 

By E. S. BiELER, M,Sc., 1851 Scholar of McGill University. 
(Communicated by Sir E. Eutherford, F.R.S.— Eeceived April 22, 1921.) 

[Plate 1.] 

The investigation of the currents induced in a conductor by the passage of 
a mass of magnetic material over it is a matter of some considerable interest. 
In the simple case where the mass is in the form of a spherical shell, the 
problem can be investigated mathematically, and the motion of the spot of a 
galvanometer connected to such a conductor can be deduced. It has occurred 
to the writer to compare the results of a mathematical investigation with those 
of experiments carried out in the laboratory on a small scale. It is the object 
of the present paper to give a summary of both the theory and the experi- 
mental confirmation. 

Theoretical Investigation. 

A. The Induced Electromotive Force, 

Consider an iron shell (fig. 1) of internal and external diameters a and h, 
and of permeability //., passing with velocity v over an infinitely long straight 
conductor at depth d below the level of its centre. Take axes through the 
centre, the Y-axis being vertical and the X- and Z-axes being horizontal, the 
one perpendicular, the other parallel to the direction of the conductor. Let 
/, m, n be the direction cosines of F, the earth's field, and let a be the angle 
between the Z-axis and the direction of motion of the shell. 

By the application of spherical harmonics it is easily proved that the 

magnetic potential at a point {x, y, z) at a distance r from the origin is 

given by 

V = 'F{lx-{-my'\-nz) (1 + X/r^), 

where \ = -(&^^a^) ~ 1?,^"^ Vxfer"^\/^.x' 

The term F {Ix + my + nz) is due to the earth's undisturbed field. The other 
term gives the potential of the disturbance. 

Instead of considering the conductor as stationary and the shell as moving- 
past it, it is simpler to reduce the shell to rest, and to consider the conductor 
as cutting the field due to the disturbance, that is, the field derived from the 

potential 

^f __ T^^ {Ix + my + nz) 
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At time t let x be the horizontal distance of the shell from the conductor. 
In the interval St, x has increased to x + Ba;, and the number N" of tubes of 
induction cut by the conductor in the interval is given by 
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The E.M.F. induced in the conductor is therefore given by 

= —dl^/dx , dx/dt 



■00 



3V' 
dy 



^ dz] V Hin oc 



FXz^sin a 



LM"- 



Sy(lx-{- my + nz)' 



/vti 



dz 



c.^^ . m(o(P—d^)'-^2lxd 



(1) 



E 2 



52 Mr. E. S. Bieler. On the Currents induced in a 

B. The Deflection produced in a critically damped Galvanometer, 

It will be assumed that the galvanometer is critically damped, a condition 
which can always be secured in a moving coil instrument by introducing a 
suitable shunt or series resistance. 

Let A be the total area of the turns of a galvanometer coil suspended in a 
field of strength H. Let E be the total resistance in the galvanometer 
circuit, K the moment of inertia of the suspended system, and r the 
directing couple of the suspension for unit angular deflection. 

If 9 represent the angular deflection of the moving system, i the current at 
any time, and E the impressed electromotive force, the two following equa- 
tions are obtained at once from the mechanical and electrical conditions 

respectively, 

¥.dWldt^'-KRi-\-Te = 0, (2) 

E = Ei + AH^(9/(^jJ. (3) 

Eliminating i between (2) and (3), we obtain 

dt^ EK dt'^'K EK ' ^ ^ 

The galvanometer will be critically damped, provided that E has been so 

adjusted that 

(A^HVEKf =:4t/K, 

that is, if E = A2I^V2v'(tK). 

If we now write n ^ A^H^/^EK = ^(t/K), equation (4) reduces to 

dW/dt^ + 2 ^ de/dt 4- nW = (2 n/AK) E, 
or {I)+nfe = i2n/Ali)'E, (5) 

The values of the coefficients n and 2^/ AH are easily obtained by two 
simple observations on the galvanometer. 

For a steady deflection, (5) reduces to nW = (2^/ AH) E, that is, 

e/B = 2/nAli, (6) 

Thus, 2/?^AH is the sensitivity in radians per unit E.M.F. for a steady 

deflection. 

Also, if the galvanometer is allowed to swing on open circuit, the motion is 

given by 

d^e/dt^ + nW = 0, 

and consequently the period T by 

T = 27r/n. (7) 

Thus, by means of (6) and (7), we can evaluate the constants of (5) for the 
voltage sensitivity and the period of the galvanometer. 

If we can solve (5) in the case where E is given by equation (1), we are 
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able to deduce the actual motion of the galvanometer spot due to the electro- 
motive force induced in a conductor by the passage of a spherical shell of 
iron over it. 

Substituting in (5) for the value of E given by (1), we obtain 

(D + ?i)2 6^ = - (4 7i/ AH) Y\v sin u (m {x^ - ^^) - 2 lnd)l{x^ 4- dy, 
which, on replacing x by its value v sin a^ = v't, reduces to 

{D-^nfe = ^ {4:71 1 AH) FXt;' (m {v'H^ - cP) - 2 Wtd)l{v'V? + d^ 

= -^M(m(^2_^2)_2%)/(^2^/)2, (8) 

where j^ = djv' and M = 47iFX,/AH'y'. 

Now, it is easily found by the method of " Variation of Parameters,"^ that 
the solution of 

is (9 = {t\e''^f{t)dt'-'lte''*f{t)dt} er'^K (9) 

Since /(?5) is of the form 

fit) = -M(m(^2_^2)_2%)/(?{2+i?7, 
the solution will contain the four unknown integrals 

^3 { t^ ( t f 1 

It would be no simple task to evaluate these for various values of n and p 
over a considerable range of t, and we shall have to restrict ourselves to the 
special case where p = Ijn. 

If we write nt = y, the four above integrals reduce to 



I. = 



is = 



f^ y^ 

J-^^^^""^' 



y y 

'^ 0} fly 



I4 



where the limits have been put in to satisfy the conditions ^ = 0, ^ = at 
infinity. 

The evaluation of these integrals will be made somewhat simpler if we 
define two other integrals, 

■* Forsyth, ' Differential Equations,' 4th ed., § 65. 
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The first three integrals may be written down in terms of the simpler 
integrals I4, Is, and Te, 



3 



1 WT^f^'^'y = -* L """^l) = ^''-2W^y 



T T T 

±2 — Iq — 14, 

T- -r ~r 111 ^ 

Substituting 3/ for nt in (9) and (10) and writings = 1/% we obtain 

= {(y/n')leycl>(i/)dy^(l/n')lyey<l>(y)dy}e~^J, (11) 

where ^(^) = --.^^2M!?i(^^l^I^ 

If we use this value of (j> {y) in (11), we finally get for 9 in terms of I4, T5, 
and Tg, 

6 = M(m{[23/l4--(s^+l)l6 + l5]6-^ + l/(^^ + l)} 

+ U[(^-2)Ie + 2Ij6~^--2//(|/^ + l)}) 
= m6>i + %. (12) 

The values of I4, I5, and Ig at intervals of 0*6 from —4*8 to 7*2 are given in 
Table I and were calculated as follows : — The integrals between the limits 
— 00 and were first worked out numerically from formulae involving the 
sine- and cosine-integrals, given by Bierens de Haan.* The integrand was 
then computed for values of the variable at intervals of 0'15 and the integrals 
evaluated at intervals of 0*6 by the use of Simpson's rule for five equidistant 
ordinates^ viz., 

ydx= ^;?,(7^i + 323/'2+123/3 + 32y4+7^5), 



1: 



-2/1 



where h is the interval between successive ordinates, and yi, y^, yz, y4. and y^ 
are the values of the integrand at each ordinate. The computations were 
carried out with the aid of an excellent Brunsviga calculating machine 
kindly lent by Mr, Eichmond, of the Ballistic Office. 

In addition to the integrals I4, I5 and Iq, Table I gives the values of OijM. 

* * Yerhand. der Kon. Akad. van Wet.,' Part lY, 1858, Table 130, formulse 4 and 5 ; 
Table 131, formula 9 ; and * Nouvelles Tables d'lnt^grales Definies,^ 1867, Table 91, 
formulse 7 and 8 ; Table 92, formula 6. 

(Formula 9, Table 131, of the * Yerhandelingen ' is correct as it stands, namely, 



'°° Q-PX 



—T dx — 1/29^ [ci {pq) . sin (p^) ~ si {pq) . cos {pq) + ^Tr . cos (pq) 



- pq {ci (pq) . cos (pq) + si (pq) , sin (pq) ~ ^w . sin (pq)}} 

The minus sign before the pq has been changed by mistake to a plus sign in the 
'' Nouvelles Tables," Table 92, formula 6.) 
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and O2/M.. These deflections are also plotted in the graphs of figs. 2 and 3. 
In each case the dotted curve represents the E.M.r., and the full curve the 
resulting galvanometer deflection. Fig. 2a corresponds to the case Z = 0, 
m = 1 ; fig. 2& to Z == 1, m == ; and fig. 3 to / = sin S/i/(2), m = cos S, 
where S is the dip, here equal to 75°. 

Table I. 
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Experimental Oonfirmation. 

The truth of the above results was tested by experiments carried out on a 
small scale in the Electrical Measurements Laboratory at McGill University 
by the use of a 16-lb. spherical iron shot and a coil of many turns connected 
to a reflecting galvanometer, whose deflections were registered photo- 
graphically. 

The coil, on a square frame of 125 cm. side, had 200 turns of No. 28B. and 
S. copper wire, wound close together, to give a circular section of about 
1 cm. diameter. Its resistance, 221*5 ohms, was made up to 224 ohms, the 
critical damping resistance of the galvanometer. 

The galvanometer was a moving coil instrument of 8 seconds free period, 
whose sensitivity, when in series with its critical damping resistance, was 
1*63 mm. per microvolt at 1 metre. 

The iron sphere, an iron shot kindly lent by Spalding and Co., proved to 
have excellent magnetic qualities. It exhibited practically no permanent 
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magnetism after lying in the same position in the earth's field for several 
days. 

The galvanometer deflection was recorded by means of a camera consisting 
of a box with a long horizontal slit. Metal cylinders, actuated by clockwork, 
unwound a roll of sensitised paper past the slit. A time-mark was made 
every second by a small pea-lamp connected to a clock. In addition, the 
pea-lamp circuit was made for an instant as the shot came past three 
contacts, one directly opposite the side of the coil, and the other two 40 cm. 
before and after, respectively. 

Experiments were carried out for three different positions of the coil. 

A. Vertical Comijonent. — The shot was made to travel horizontally in a 
magnetic east and west direction over one side of the coil, pointing magnetic 
north and south. The horizontal component of the earth's field acted along 
the side of the coil, and the vertical component alone was effective. In 
this case, / = and m = sin 8 = sin 75^ in formulae (1) and (8) (see fig. 2a). 




Since the theoretical investigation had only been carried out for the case 
p = l/?i, it was necessary to make sure that the velocity of the shot past the 
coil was such as to secure this condition. The parameters p and n are 
connected with the period of the galvanometer, the height of the shot above 
the coil, and the velocity at right angles to the coil, by the two relations 
T = 27r/?i (equation 7) and p = cljv' (equation 8). 
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Eliminating p and n between these three equations, we have the relation 

v' = 2ml/T, (13) 

which gives the required velocity. 

Plate la gives a record of the galvanometer deflection in this case for 
d = 14'6 cm. The thick black line is the deflection record, the thin grey 
lines on either side being only ghosts. The deflection agrees very well with 
the curve of fig. 2a, although the shot velocity was a little low, 10*2 cm. per 
second instead of 11-6. The maximum deflection measured on the record is 
17*0 mm., which is in fair agreement with 14'3 mm., the deflection calculated 
from formula (12). The difference is probably due to slight permanent 
magnetisation of the shell. 

With the coil in a horizontal position, trouble was experienced at first, due 
to large disturbances presumably caused by induction from the tramways in 
the neighbourhood. These disturbances were completely cut down by con- 
necting in series with the coil, and in opposition to it, three turns of wire 
laid round the laboratory, and enclosing a total area equal to that of all the 
turns in the coil. 

B. Horizontal Component. — In this case the coil was placed in a vertical 
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plane pointing magnetic north and south, and the shot was moved up 
and down past one of its horizontal sides. The horizontal component 
of the earth's field was directed along the side of the coil, and the 
vertical component, acting along the direction of motion of the shot, 
played the part of a virtual horizontal component. This corresponds to 
I = sin S = sin 75° and m = in formulae (1) and (8) (see fig. 2&). 

Here, again, the velocity was determined by means of formula (13). JSTo 
noticeable tramway disturbances were observed in this position of the coil. 
Plate lb gives a record of the galvanometer deflection for d = 8*2 cm. 
There is good agreement with the calculated deflection shown in fig. 2&, 
although the velocity of the shot was somewhat low. The maximum 
deflection is 35 mm., against 32*2 mm. calculated from formula (12). 

C. Both Components Comhined. — The coil was again placed in a vertical 
position, but pointing in a magnetic east and west direction, with its sides 
inclined at 45° to the horizontal. Here, Z = sinS/ \/(2) = sin 75°/ \/(2), 
m = cos B. As may be seen from equation (1), the E.M.F. due to the com- 
ponent parallel to the direction of motion of the shot changes sign with x, 
whereas that due to the component perpendicular to both the side of the coil 
and the direction of motion does not. The form of the E.M.F. curve and of 
the deflection curve will thus depend on the direction of motion of the shot. 
This is illustrated by the curves of fig. 3 and the record of Plate Ic. In 
fig. 3a, and the first deflection in the plate, the positive maximum comes just 
before the zero-line, and is approximately equal to the negative maxinmm. 
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In fig. 35, and the second deflection in the plate, the positive maximum 
comes just before the zero-line, and is over twice the negative maximum. 




Conclusion. 

We have now shown that the theory agrees with experimental evidence, 
and we may draw further conclusions as to the way in which the galvanometer 
deflection varies with different factors. 

The galvanometer deflection is proportional to a factor, M, given by 



ivi — — — -_ z . 



2i |7, ^ 1 Jtt 



Thus, if the velocity is that given by D'-27rd/T (equation 13), the 
deflecfcion will be proportional to — 

(a) The voltage sensitivity, 2/AH%, of the galvanometer; 

(h) The earth's field, F ; and 

(c) The constant X, which depends on the amount and distribution of 
magnetic material in the shell, 
and will be inversely proportional to — 

(a) The depth d of the conductor below the level of the shell ; and 

(&) The period T of the galvanometer. 

Further, it has been found experimentally, and pointed out in the case of 
Plate 1, a and &, that a small variation in the velocity produces only a small 
effect on the deflection. The above results will therefore hold approximately 
even when the velocity is not that given by equation (13). 

The theory also throws light upon the connection that may be expected 
between the deflection of a critically damped galvanometer, and the E.M.F. 
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variation which caused it. Two points appear at once on examination of the 
theoretical curves of figs. 2 and 3 : — 

{a) The galvanometer deflection does not give a perfectly faithful, but a 
slightly distorted, reproduction of the E.M.F. curve. 

(6) The galvanometer deflection lags behind the E.M.F. by an amount here 
equal to about one-fourth of the period. 

Both these deviations from a true representation of the E.M.F. variation 
are due to the effect of the inertia of the moving system of the galvanometer. 



On some Neiv FormidcB for the Numerical Calculation of the 

Mutual Induction of Coaxial Circles, 

By Louis V. King, D.Sc, Macdonald Professor of Physics, McGill 

University, Montreal. 

(Communicated by Prof. A. S. Eve, F.R.S. Received April 13, 1921.) 

Section 1. — Note on the Direct Nmnerical Oalculation of Elliptic Integrals. 

The importance of realising rapid and accurate methods of calculating the 
elliptic integrals, now denoted by 

F (<^, h) = {^d^/A (^, k), E (<^, ^) = f A (<jb, k) d<j,, (1) 

Jo Jo 

where A (^, k) = y^ (1 — F sin^ ^), 

was first remarked by Euler* (1766), although it was not until several years 
later that Landenf (1775) discovered in geometrical form the transformation 
which forms the basis of existing methods of numerical calculation of the 
elliptic integrals, 

A method of successive transformations for the ultimate reduction of the 
algebraic forms of the elliptic integrals to elementary forms was published 
by Lagrange:!: in 1784-5. This memoir contains an exposition of the scales 
of arithmetico-geometrical means and discusses their use in calculating the 

*^ Enler, 'Novi Coinm. Acad. Sc. Petrop.,' vol. 10 (1766). 

t Landen, 'Phil. Trans. Roy. Soc.,' vol. 65, p. 283 (1775); also, 'Mathematical 
Memoirs,' London, 1780. 

J Lagrange, " Sur une nouvelle m6thode de calcul integral pour les differentielles 
afFectees d'un radical carr6 sous lequel la variable ne passe pas le quatrieme degr6," 
'Mem. de FAcad. Koy. des Sc. de Turin,' vol. 2 (1784-5) ; ' CEuvres,' vol. 2, pp. 253-312 
(Gauthier-Villars, Paris, 1868). 



